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Abstract. We numerically investigate the rheological properties of jammed frictionless granular materials under an oscillatory 
shear in terms of our simulation of the distinct element method. It is demonstrated that the constitutive relation between the 
shear stress and the shear strain of grains strongly depends on the amplitude and the frequency of the oscillatory shear. For 
a small amplitude region, the granular material behaves as a visco-elastic material characterized by the Kelvin- Voigt model, 
while it behaves as a yield stress liquid for a large amplitude region. We find that the rheology of grains under the large 
amplitude oscillatory shear can be described by a phenomenological constitutive model. 
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INTRODUCTION 

Athermal disordered materials such as colloidal suspen- 
sions, foams, and granular materials behave as dense liq- 
uids when the density is lower than a critical value, while 
they behave as amorphous solids when the density ex- 
ceeds the critical value [1, 2, 3]. This rigidity transition 
is known as the jamming transition. 

The jamming transition has attracted many physicists 
since Liu and Nagel indicated its similarity to the glass 
transition [4]. Recently, the rheology of jammed materi- 
als under steady shear has been extensively studied. It is 
remarkable that there exist critical scaling laws for the 
rheological transition of frictionless grains for the pres- 
sure, the elastic modulus, and the shear stress, similar to 
those in continuous phase transitions, though the number 
of contact points is discontinuously changed at the jam- 
ming point [5, 6, 7, 8, 9, 10, 1 1, 12, 13, 14, 15, 16, 17, 18, 
19, 20, 21]. This mixed phase transition becomes the dis- 
continuous transition associated with the hysteresis loop 
for frictional granular materials [22, 23, 24]. 

On the other hand, there exist only a few studies on 
the rheological properties of jammed materials near the 
jamming transition point under an oscillatory shear. Re- 
cently, a numerical simulation of grains in a connected 
spring network suggested that there is a critical scaling 
law for the complex shear modulus under an oscillatory 
shear [25]. However, the model in Ref. [25] underesti- 
mates roles of rattlers or the change of connecting net- 
works, and thus, it may be unsuitable to describe real- 
istic jammed materials. Moreover, as long as we know, 
nobody has systematically analyzed rheology of granu- 
lar materials under the oscillatory shear. 



In this paper, we numerically investigate the rheolog- 
ical properties of jammed frictionless granular materials 
under the oscillatory shear. In the next section, we ex- 
plain our setup and models. We demonstrate that the rela- 
tion between the shear stress and the shear strain strongly 
depends on the amplitude and the frequency of the oscil- 
latory shear, and the rheological properties are consis- 
tent with predictions by a phenomenological constitutive 
model in the third section. Finally, we discuss and con- 
clude our results in the last section. 



MODEL AND SETUP 

Let us simulate a two-dimensional frictionless granular 
assembly in a square box with side length L in terms of 
the distinct element method (DEM). The system includes 
N grains, each having an identical mass m. The position 
and velocity of a grain i are respectively denoted by r, 
and V/. Our system consists of grains having the diam- 
eters 0.7<7o, 0.8(7(), 0.9(7o, and Co, where the number of 
each species of grains is N/4. The contact force between 
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grains ; and j consists of the elastic part fj. and the dis- 
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sipative part ft. , which are respectively given by 
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where ny is ny = ry/|ry| with the normal elastic con- 
stant k, the viscous constant 77, the diameter (7, of grain i, 

ry = r ; -Tj, (Ty = (C7; + Oj)/2, and Vy = (v; - Vj) ■ ny. 

Here, ®(x) is the Heaviside step function characterized 
by &(x) = 1 for x > and 0(jc) =0 for otherwise. 



In this paper, we apply an oscillatory shear along the 
y direction under the Lees-Edwards boundary condition 
[26]. As a result, there exists macroscopic displacement 
only along the x. The time evolution of such a system, 
known as the SLLOD system, is given by 

~j~ = - + r(t)y.^ (3) 

at m 

with the peculiar momentum p,- and the unit vector par- 
allel to the x-direction e x . The shear rate and the shear 
strain are, respectively, given by 

f(t) = fofl)sin(a»), (5) 
7(0 = ro{l-cos(a)f)} (6) 

with the amplitude 70 and the angular frequency CO. In 
this paper, we investigate the shear stress 5: 

(7) 

where (•) represents the ensemble average. 

In our simulation m, Go, and k are set to be unity, 
and all quantities are converted to dim ensionless forms, 
where the unit of time scale is y/m/k. We use the vis- 
cous constant 77 = 1 .0 This situation corresponds to the 
constant restitution coefficient e = 0.043. We use the 
leapfrog algorithm, which is second-order accurate in 
time, by using the time interval At = 0.2. The number 
N of the particles is 4000. We fix the volume fraction 
of the system = 0.67, which is larger than the critical 
fraction <pj = 0.645 under a steady shear [19]. 
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FIGURE 1. Shear stress 5(f) as a function of shear strain y(t) 
for the amplitude y = 1(T 4 with (O = KT 1 , 1(T 2 , and 1(T 3 . 
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FIGURE 2. Shear stress 5(f) as a function of shear strain 
y(f) for the amplitude 70 = 10° with a = 1CT 1 , 1(T 2 , 1(T 3 , 
and 1(T 4 . 



RESULTS 

Figure 1 shows the shear stress 5(f) as a function of 
the shear strain 7(f) for the amplitude 70 = 10~ 4 with 
the frequency CO = 10~', 10~ 2 , and 10~ 3 . For co = 10~\ 
the trajectory draws an ellipse. The width of it decreases 
as co decreases. Finally, the trajectory becomes a single 
straight line for small co. This is a typical property of 
the visco-elastic materials described by the Kelvin- Voigt 
model, whose constitutive equation is given by 



5 = Gy+ 777 



(8) 



with the shear modulus G and the shear viscosity 77 . The 
application of the Kelvin- Voigt model might be natural 
because the shear amplitude is so small that the change 
of the contact network is negligible. 

On the other hand, the rheology of jammed grains dif- 
fers from that of the Kelvin- Voigt model for the larger 



amplitude 70 = 1.0. In Fig. 2, we plot the shear stress 
5(f) as a function of the shear strain 7(f) for the ampli- 
tude 70 = 10° with CO = lO-SlO^VO , and 10~ 4 As 
the frequency decreases, the width of the trajectory de- 
creases, but it still remains even in the smallest CO. For the 
smallest frequency CO = 10~ 4 , the shear stress 5 shows a 
linear elastic dependence on the shear strain in the region 
< 7 < 0.2 with 5 > and 1 .8 < 7 < 270 with 5 < 0, 
while |5| is almost constant for the other region. The ex- 
istence of the width in the limit co —> is a typical prop- 
erty of yield stress fluids [27]. 

In order to describe the rheological properties of yield 
stress fluids, such as aqueous foams [28], the elastic 
Bingham model has been used [27, 29]. The shear stress- 
strain relationships for the elastic Bingham model are 
given by 



5 = 



GJe, 
Gy c - 



\Ye\ < 7c 
T77, bfel > Yc 



(9) 



where y c is the yield strain and Je is the recoverable elas- 
tic strain Je- The elastic strain Ye is obtained by inte- 
grating the shear rate with respect to time, but saturates 
at \Ye\ = Yc- When the shear stops and the direction of 
it is reversed, the accumulated elastic strain is recovered 
and the response is linear elastic until the material is re- 
yielded again for \y E \ = y c - 

In the elastic Bingham model, the shear stress under 
steady shear with the shear rate y is given by 



0.04 



S = Gy c + r\y, 



(10) 



but this is inconsistent with the constitutive equation of 
frictionless granular materials under a steady shear. In- 
deed, the shear stress S of the jammed frictionless gran- 
ular materials under the steady shear satisfies a critical 
scaling law given by [13, 14] 



s=(<i>-(i>j)y^F[y/(<j>-<j> j )y^yy 

where the scaling function F(x) obeys 
limF(x) = const., 



lim F(x) 
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The critical exponents y$ and y y for the interaction force 
in Eq. (1) is approximately given by 



» = 1, y r = 2/5, 



(14) 



which are also phenomenologically derived [13, 14]. 
Hence, the shear stress S of the jammed frictionless gran- 
ular materials under steady shear is expected to satisfy 

S = S c + ri'f/ 5 (15) 

with a constant rj' and the yield stress S c , which differs 
from that of the Kelvin- Voigt model (10). 

In order to characterize the rheology of grains as in 
Fig. 2, we propose a new phenomenological constitutive 
equation, where the shear stress S is given by 

Gy E + r]'f 15 , \Ye\<Yc 



S = 



Gy+ri'y 



'y 2 /5 



\Ye\ > Yc- 



(16) 



Here, we assume that the shear stress consists of the 
dynamic part given by r\'y 2 l s and the static part, which 
is proportional to the elastic strain y E until the strain 
exceeds the yield strain y c , where the change of the 
contact network appears. 

In Fig. 3, we plot the shear stress S(t) as a func- 
tion of the shear strain y(t) obtained from the phe- 
nomenological model (16) for the amplitude 70 = 10° 
with co = lO^lO^lO" 3 , and 1(T 4 . Here, we choose 
the parameters as G = 0.02, r/' = 0.06, and y c = 0.06. 
From the comparison of Fig. 2 with Fig. 3, we find that 
our phenomenology given by Eq. (16) provides semi- 
quantitatively accurate behavior of DEM shown in Fig. 
2. 
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FIGURE 3. Shear stress S(t) as a function of shear strain 
y(t) obtained from the phenomenological model (16) for the 
amplitude y = 10° with a) = 10 J , 10 2 , and 10~ 3 . 



DISCUSSION AND CONCLUSION 

It is known that the rheological properties of granular 
materials under a steady shear drastically depend on 
the volume fraction near the jamming transition point. 
[9, 10, 12, 13, 14, 15, 19]. We also expect that the 
rheological properties of granular materials under an 
oscillatory shear strongly depend on the volume fraction 
0. We will report such a -dependence of the constitutive 
equation elsewhere. 

In this paper, we restrict our interest to the frictionless 
particles. When the particles have friction, the rheologi- 
cal properties under oscillatory shear should be changed. 
In fact, the critical properties and the critical fraction de- 
pend on the friction coefficient for frictional granular par- 
ticles under steady shear [22, 24]. The rheological prop- 
erties in the dynamical systems of frictional particles un- 
der oscillatory shear will be discussed in our future work. 

When the amplitude 70 becomes small, the constitutive 
model given by Eq. (16) converges to 



S = Gy+r]'f l \ 



(17) 



which differs from the Kelvin- Voigt model (8). For small 
70, we expect that the rheology is characterized by the 
Kelvin- Voigt. Hence, we should improve the constitutive 
model in order to describe the rheological properties in 
wide range of 70. 

In conclusion, we have numerically studied fric- 
tionless granular materials under oscillatory shear, and 
demonstrated that the rheology depends on the ampli- 
tude of the shear. We also proposed a new phenomeno- 
logical constitutive equation, and demonstrated that the 
rheology for the large amplitude is well described by our 
phenomenological model. 
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